Ancilla Assisted Quantum State Tomography in Many-Qubit Registers 
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The standard method of Quantum State Tomography (QST) relies on the measurement of a set of 
noncommuting observables, realized in a series of independent experiments. Ancilla Assisted QST 
(AAQST) proposed by Nieuwenhuizen and co-workers (Phys. Rev. Lett., 92, 120402 (2004)) greatly 
reduces the number of independent measurements by exploiting an ancilla register in a known initial 
state. In suitable conditions AAQST allows mapping out density matrix of an input register in a 
single experiment. Here we describe methods for explicit construction of AAQST experiments in 
multi-qubit registers. We also report nuclear magnetic resonance studies on AAQST of (i) a two- 
qubit input register using a one-qubit ancilla in an isotropic liquid-state system and (ii) a three-qubit 
input register using a two-qubit ancilla register in a partially oriented system. The experimental 
results confirm the effectiveness of AAQST in such many-qubit registers. 
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I. INTRODUCTION 

Quantum computers have the potential to carry-out 
certain computational tasks with an efficiency that is 
beyond the reach of their classical counterparts [1]. In 
practice however, harnessing the computational power 
of a quantum system has been an enormously challeng- 
ing task [2] . The difficulties include imperfect control on 
the quantum dynamics and omnipresent interactions be- 
tween the quantum system and its environment leading 
to an irreversible loss of quantum coherence. In order 
to optimize the control fields and to understand the ef- 
fects of environmental noise, it is often necessary to com- 
pletely characterize the quantum state. In experimental 
quantum information studies, Quantum State Tomogra- 
phy (QST) is an important tool that is routinely used to 
characterize an instantaneous quantum state [1]. 

QST on an initial state is usually carried out to confirm 
the efficiency of initialization process. Though QST of 
the final state is usually not part of a quantum algorithm, 
it allows one to measure the fidelity of the output state. 
QSTs in intermediate stages often help experimentalists 
to tune-up the control fields better. 

QST can be performed by a series of measurements of 
noncommuting observables which together enables one 
to reconstruct the complete complex density matrix. In 
the standard method, the required number of indepen- 
dent experiments grows exponentially with the number 
of input qubits [3, 4]. Anil Kumar and co-workers have il- 
lustrated QST using a single two-dimensional NMR spec- 
trum [5] . They showed that a two-dimensional NMR cx- 
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periment consisting of a series of identical measurements 
with systematic increments in evolution time, can be used 
to quantitatively estimate all the elements of the den- 
sity matrix. Later Nieuwenhuizen and co-workers have 
shown that it is possible to reduce the number of inde- 
pendent experiments in the presence of an ancilla register 
initialized to a known state [6] . They pointed out that in 
suitable situations, it is possible to carry-out QST with a 
single measurement of a set of factorized observables. We 
refer to this method as Ancilla Assisted QST (AAQST). 
This method was experimentally illustrated by Suter and 
co-workers using a single input qubit and a single ancilla 
qubit [7] . Recently Peng and coworkers have studied the 
effectiveness of the method for qutrit-like systems using 
numerical simulations [8]. 

Single shot mapping of density matrix by AAQST 
method not only reduces the experimental time, but also 
alleviates the need to prepare the target state several 
times. Often slow variations in system Hamiltonian may 
result in systematic errors in repeating the state prepa- 
ration. Further, environmental noises lead to random 
errors in multiple preparations. These errors play impor- 
tant roles in the quality of the reconstruction of the target 
state. Therefore AAQST has the potential to provide a 
more reliable way of tomography. 

In this article we first revisit the theory of QST and 
AAQST and provide methods for explicit construction of 
the constraint matrices, which will allow extending the 
tomography procedure for large registers. An important 
feature of the method described here is that it requires 
only global rotations and short evolutions under the col- 
lective internal Hamiltonian. We also describe nuclear 
magnetic resonance (NMR) demonstrations of AAQST 
on two different types of systems: (i) a two-qubit input 



register using a one-qubit ancilla in an isotropic liquid- 
state system and (ii) a three-qubit input register using a 
two-qubit ancilla register in a partially oriented system. 
In the following section we briefly describe the theory 
of QST and AAQST. In section III we describe experi- 
mental demonstrations and finally we conclude in section 
IV. 



II. THEORY 



jth bit. Here v = v x 2 n - 2 + v 2 2 n -' i + h f„_i2° is the 

value of the n — 1 bit binary string (i/ 1; i>i, • • • , v n -\) and 
v can take a value between and 7 = N/2 — 1. The 
real and imaginary parts of an NMR signal recorded in 
a quadrature mode corresponds to the expectation val- 

n 

ues of transverse magnetization observables ^2 o~j X and 

n 

J2 a jy respectively [9]. 
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A. Quantum State Tomography 

We consider an n-qubit register formed by a system of 
n mutually interacting spin- 1/2 nuclei with distinct reso- 
nance frequencies u>i and mutual interaction frequencies 
2irJij. The Hamiltonian under weak-interaction limit 

(27rJjj < \Ui - LOj 



consists of the Zeeman part and 



spin-spin interaction part, i.e. 
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respectively, where at, and a J z are the z-components of 
Pauli operators of ith and jth qubits [9]. The set 
of N — 2™ eigenvectors {|"z 1 7n 2 • • -to„)} of the Zee- 
man Hamiltonian form a complete orthonormal compu- 
tational basis. We can order the eigenvectors based on 
the decimal value m of the binary string (mi • • • m n ), i.e., 

m = mi2 n - x H hm„2°. 

The general density matrix can be decomposed as 
t/N + ep where the identity part is known as the back- 
ground, the trace-less part p is known as the deviation 
density matrix, and the dimensionless constant e is the 
purity factor [10]. In this context, QST refers to complete 
characterization of the deviation density matrix, which 
can be expanded in terms of N 2 — 1 real unknowns: 

N-2 

p=J2p™m(\m)(m\-\N-l)(N-l\) 

m=0 

N-2 W-l 

+ E E {-Rmm'(|m)(m'| + |m'><m|) 

m=0 m'=m+l 

+ iS mml (\m)(m'\-\m')(m\)}. (2) 

Here first part consists of N — 1 diagonal unknowns 
Pmm with the last diagonal element Pn-i,n-i being con- 
strained by the trace-less condition. R and S each con- 
sisting of (N 2 — N)/2 unknowns correspond to real and 
imaginary parts of the off-diagonal elements respectively. 
Thus a total of N 2 — 1 real unknowns needs to be deter- 
mined. 

Usually an experimental technique allows a par- 
ticular set of observables to be measured directly. 
To explain the NMR case, we introduce n-bit bi- 
nary strings, j v = v x v 2 - ■ -v^iOvj ■ ■ ■ v n -i and j' v = 
V1V2 ■ ■ ■ v j-\^ v j ' ' ' v n-\ differed only by the flip of the 



The background part of the density matrix neither 
evolves under unitaries nor gives raise to any signal, and 
therefore we ignore it. Under suitable conditions (when 
all the transitions are resolved), a single spectrum di- 
rectly yields nN matrix elements {Rj vt j> , Sj vi j> } as com- 
plex intensities of spectral lines. These matrix elements 
are often referred to as single quantum elements since 
they connect eigenvectors related by the flip of a sin- 
gle qubit. We refer the single-quantum terms Rj„,j' and 
Sj v ji respectively as the real and imaginary parts of z^th 
spectral line of jth qubit. Thus a single spectrum of an 
n-qubit system in an arbitrary density matrix can yield 
nN real unknowns. 



In order to quantify the remaining elements, one relies 
on multiple experiments all starting from the same ini- 
tial state p. The fcth experiment consists of applying a 
unitary Uk to the state p, leading to p^ = UkpU* k , and 
measuring the single-quantum spectrum {RS '., , S, ,, }. 
From eqn. (2) we obtain 

R Tlti = E 4? ( m )Pmrn + 

III 

E C jv ( TO ' m ') R ™rn> + e)J (m, m')S mrn > , 

m,Tn'>7n 

S jt,Jl = E b fJ( m )Pmrn + 
m 

^ djj (m, m!)R mm , + fjjj (to, m')S mm i ,(3) 

in terms of the unknowns p mm ' an d the known real con- 
stants {a, • ■ • , /}: 

afj (m, to) + ibfj (to, m) = (jv\U k \m}(m\Ul\j' v ) - 

<ju\U k \N-l){N-l\Vt\? v ), 

cf}(m,m') + idfj(m,m') = {j„\U k \m)(m'\Ul\j' v ) + 

(j v \U k \m')(m\Ut\j' u ), 

efj(m,m') + if£ ) (m,m') = i{j v \Uk\m)(m'\u£\Jl) - 

i(ju\U k \m')(m\Ut\Jl) (4) 

[11]. After K experiments, we can setup the matrix equa- 
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Here the left column vector is formed by the TV 2 — 1 
unknowns of p: diagonal elements in the top, real off- 
diagonals in the middle, and imaginary off-diagonals in 
the bottom. The right column vector is formed by KnN 
numbers - the real and imaginary parts of the experi- 
mentally obtained spectral intensities ordered according 
to the value of the binary string v, the qubit number j, 
and the experiment number k. The KnN x (TV 2 — f) 
dimensional constraint matrix is of the form 

M = 
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Note that each column of the constraint matrix corre- 
sponds to contribution of a particular unknown element 
of p to the various spectral intensities. 



By choosing the unitaries {£/&} such that rank(M) > 
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FIG. 1. Minimum number of independent experiments re- 
quired for QST (with zero ancilla) and AAQST. 



TV 2 — 1 (the number of unknowns), eqn. (5) can be solved 
either by singular value decomposition or by Gaussian 
elimination method [11]. Fig. 1 illustrates the minimum 
number (K) of experiments required for QST. As antic- 
ipated, K increases rapidly as 0(N/n) with the number 
of input qubits. In the following we describe how it is 
possible to speed-up QST, in the presence of an ancilla 
register, with fewer experiments. 



B. Ancilla Assisted QST (AAQST) 

Suppose the input register of n-qubits is associated 
with an ancilla register consisting of n qubits. The di- 
mension of the combined system of h = n + h qubits is 
N = NN, where N = 2 n . For simplicity we assume that 
each qubit interacts sufficiently with all other qubits so 
as to obtain a completely resolved spectrum yielding hN 
real parameters. Following method is applicable even if 
there are spectral overlaps, albeit with lower efficiency 
(i.e., with higher number (K) of minimum experiments). 
Further for simplicity, we assume that the ancilla regis- 
ter begins with the maximally mixed initial state, with 
no contribution to the spectral lines from it. Otherwise, 
we need to add the contribution of the ancilla to the final 
spectrum and the eqn. (5) will become inhomogeneous. 
As explained later in the experimental section, initializa- 
tion of maximally mixed state can be achieved with high 
precision. Thus the deviation density matrix of the com- 
bined system is p = p <Z> 1/N. Now applying only local 
unitaries neither leads to ancilla coherences nor trans- 
fers any of the input coherences to ancilla. Therefore 
we consider applying a non-local unitary exploiting the 
input-ancilla interaction, 



N-l 



Uk = Vj2 Uka®\a)(a\ 



(7) 



a=Q 



where Uka is the kth unitary on the input register depen- 
dent on the ancilla state \a) and V is the local unitary 
on the ancilla. The combined state evolves to 



~ p W 



UkP&t 



4 E Pmm'U ka \m)(m'\Ul®V\a)(a\Vt. (8) 
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We now record the spectrum of the combined system cor- 

n 

responding to the observable Yl &j X +i<jj y . Each spectral 

line can again be expressed in terms of the unknown el- 
ements of the ancilla matrix in the form given in eqn. 
(3). The spectrum of the combined system yields hN 
linear equations. The minimum number of independent 
experiments needed is now 0(N 2 /(hN)). Since we can 
choose N > N, AAQST needs fewer than 0(N/ri) ex- 
periments required in the standard QST. In particular, 
when hN > N 2 . a single optimized unitary suffices for 
QST. Fig. 1 illustrates the minimum number (K) of ex- 
periments required for various sizes of input and ancilla 
registers. As illustrated, QST can be achieved with only 
one experiment, if an ancilla of sufficient size is provided 
along with. 



C. Building the constraint matrix 



The major numerical procedure in AAQST is obtain- 
ing the constraint matrix M. For calculating the con- 
straint coefficients c r , , one may utilize an elaborate de- 
composition of C/fc using numerical or analytical methods. 
Alternatively, as described below, we can use a simple al- 
gorithmic approach to construct the constraint matrix. 

First imagine a diagonal state p for the ancilla register 
(eqn. (2)) with p o = 1 an d Pmm = for all other 1 < 
m < A -2, Rmm> = S mm i = 0. Applying the unitary U k 
on the composite deviation density matrix 5 = p® 1/N, 
we obtain all the spectral intensities (using eqn. (3)) 



? (fe) 



a*, (0,0) = R)li JV „ 6^(0,0) = 5. 



?(fe) 



(9) 



Thus the spectral lines indicate the contributions only 
from poo (and pn-i.n-i)- Repeating the process with 
all the unitaries {Uk} yields the first column in M ma- 
trix (eqn. (6)) corresponding to the unknown poo- Same 
procedure can be used for all the diagonal elements p m m 
with < m < N - 2. 

To determine M matrix column corresponding to a 
real off-diagonal unknown R mm ' , we start with an input- 
register density matrix R mm ' = 1 and all other elements 
set to zero. Again by applying the unitary Uk on the 
composite density matrix, and using eqn. (3) we obtain 

c% (m, m') = R% v , , d) v (m, m') = S^ jv , . (10) 

Repeating the process with all unitaries {Uk} determines 
the column of M corresponding to the unknown R mm ' ■ 



To determine M matrix column corresponding to an 
imaginary off-diagonal unknown S mm i, we set S mm i = 
1 and all other elements to zero, and apply Uk on the 
composite state to obtain 



4,(m,m')=.R^ v , /^(m,m') = ^ v . (11) 

Proceeding this way, by selectively setting the unknowns 
one by one, the complete constraint matrix can be built 
easily. 



D. Optimization of Unitaries 

Solving the matrix equation (5) requires that 
rank(M) > A^ 2 — 1, the number of unknowns. But hav- 
ing the correct rank is not sufficient. The matrix M must 
be well conditioned in order to ensure that small errors 
in the observed intensities {RjJ- v , , SjJ ' v } do not con- 
tribute to large errors in the values of the elements p mm ' ■ 
The quality of the constraint matrix can be measured by 
a scalar quantity called condition number G(M) defined 
as the ratio of the largest singular value of M to the small- 
est [12]. Smaller the value of C(M), better conditioned is 
the constraint matrix M for solving the unknowns. Thus 
the condition number provides a convenient scalar quan- 
tity to optimize the set {Uk} of unitaries to be selected for 
QST. As explained in the experimental section, we used 
a simple unitary model U± (ti , T2) as an initial guess and 
used genetic algorithm to minimize the condition number 
and optimize the parameters (ri,T2). 

The necessary number (K) of independent experiments 
is decided by the rank of the constraint matrix and 
the desired precision. The rank condition requires that 
KnN > A 2 — 1. Introducing additional experiments 
renders the problem over-determined, thus reducing the 
condition number and increasing the precision. In the 
following section we describe the experimental results of 
AAQST for registers with (i) n = 2,n = 1, h = 3 and (ii) 
n = 3, h = 2, h = 5 respectively. 



III. EXPERIMENTS 

We report experimental demonstrations of AAQST on 
two spin-systems of different sizes and environments. In 
each case, we have chosen two density matrices for tomog- 
raphy. All the experiments described below arc carried 
out on a Bruker 500 MHz spectrometer at an ambient 
temperature of 300 K using high-resolution nuclear mag- 
netic resonance techniques. 



A. Two-qubit input, One-qubit ancilla 

Here we use three spin-1/2 19 F nuclei of iodotrifluo- 
roethylene (C2F3I) dissolved in acctone-Dg as a 3-qubit 






D D 


|Q 


11860.8 69.9 


47.4 


tl ° 


-128.3 


< 


-17379.1 1 



F1 



F2 



F3 



FIG. 2. Molecular structure of iodotrifluoroethylene, and the 
table of Hamiltonian parameters in Hz: chemical shifts (di- 
agonal elements) and J-coupling constants (off-diagonal ele- 
ments). 



system. The molecular structure and the Hamiltonian 
parameters are shown in Fig. 2. As can be seen in Fig. 3, 
all the 12 transitions of this system are clearly resolved. 
We have chosen F\ as the ancilla qubit and F 2 and F 3 
as the input qubits. QST was performed for two different 

density matrices (i) thermal equilibrium state, i.e., p\ 

i 

2 



a 2 z + erf), and (ii) state after a (7r/4)_/4 pulse applied 



to the thermal equilibrium state, i.e., p 2 



{°l 



In both the cases, 



thc first 



qubit was initialized into a maximally mixed state by 
applying a selective { / n/2) y pulse on F\ and followed by 
a strong pulscd-ficld-gradicnt (PFG) in the ^-direction. 
The selective pulse was realized by GRAPE technique 
[131. 



AAQST of each of the above density matrices required 
just one unitary evolution followed by the measurement 
of complex NMR signal. We modelled the AAQST uni- 
tary as follows: U\ = (f ) U^fa) (f ) x Umt(n), where 
Uint{r) = exp(— iHt) is the unitary operator for evolu- 
tion under the internal Hamiltonian % (see eqn. (1)) for a 
time r, and (?) rotations are realized by non selective ra- 
dio frequency pulses applied to all the spins along the di- 
rections indicated by the subscripts. The constraint ma- 
trix M had 15 columns corresponding to the unknowns 
and 24 rows corresponding to the real and imaginary 
parts of the 12 spectral lines. Only the durations {ri, t 2 } 
needed to be optimized to minimize the condition num- 
ber C(M). We used a genetic algorithm for the optimiza- 
tion and obtained C(M) = 17.3 for n = 6.7783 ms and 
t 2 = 8.0182 ms. The real and imaginary parts of the sin- 
gle shot experimental AAQST spectrum, along with the 
reference spectrum, are shown in the top part of Fig. 3. 
The intensities {-RL ',■„/ , SjJj V < } were obtained by simple 
curve-fit routines, and the matrix eqn. (5) was solved 
to obtain all the unknowns. The reconstructed density 
matrices along with the theoretically expected ones are 
shown below the spectra in Fig. 3. The fidelities of exper- 
imental states with the theoretically expected states (pi 
and p 2 ) are respectively 0.998 and 0.990. The high fideli- 
ties indicated successful AAQST of the prepared states. 
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FIG. 3. AAQST results for thermal equilibrium state pi 
(left column), and that of state p2 (right column), described 
in the text. The reference spectra is in the top trace. The 
spectra corresponding to the real part (R . . , , middle trace) 

and the imaginary part (S J . , , bottom trace) of the 19 F 
signal are obtained in a single shot AAQST experiment. The 
bar plots correspond to theoretically expected states (top row) 
and those obtained from AAQST experiments (bottom row). 
Fidelities of the states are 0.997 and 0.99 respectively for the 
two density matrices. 



B. Three-qubit input, Two-qubit ancilla 

We use three 19 F nuclei and two : H nuclei of 1-bromo- 
2,4,5-trifluorobenzene partially oriented in a liquid crys- 
tal namely, N-(4-methoxybenzaldehyde)-4- butylanilinc 
(MBBA). Due to the partial orientational order, the di- 
rect spin-spin interaction (dipolar interaction) does not 
get fully averaged out, but gets scaled down by the order 
parameter [14]. The chemical shifts and the strengths of 
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FIG. 4. Molecular structure of l-bromo-2,4,5- 

trifluorobenzene, and the table of Hamiltonian parameters 
in Hz: chemical shifts (diagonal elements) and effective 
coupling constants (J+2D) (off-diagonal elements). 
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FIG. 5. AAQST results for thermal equilibrium state, i.e., 
(p\ + a\ + (Tz)/2. The reference spectrum is in the top trace. 
The spectra corresponding to the real part (R J ■ , , middle 
trace) and the imaginary part (R-J v ,, bottom trace) of the 
19 F signal are obtained in a single shot AAQST experiment. 
The bar plots correspond to theoretically expected states (top 
row) and those obtained from AAQST experiments (bottom 
row). Fidelity of the AAQST state is 0.98. 



the effective couplings are shown in Fig. 4. As is evi- 
dent, the partially oriented system can display stronger 
and longer-range coupling network leading to a larger 
register. Here we choose the three 19 F nuclei forming 
the input register and two 1 H nuclei forming the ancilla 
register. The Hamiltonian for the heteronuclear dipolar 
interaction (between 1 H and 19 F) has an identical form as 
that of J- interaction [14]. The homonuclear dipolar cou- 
plings (among 19 F, as well as among 1 H nuclei) were small 
compared to their chemical shift differences enabling us 
to approximate the Hamiltonian in the form of eqn. (1). 



The partially oriented spin-system yields all the 80 
transitions sufficiently resolved. Again we use just one 
experiment for the complete AAQST of the 3-qubit input 
register. We modelled the AAQST unitary in a similar 
way as before: U 1 = (f ) U int (T 2 ) (f ) U int (n) where 




FIG. 6. AAQST results for the state p2 described in the text. 
The reference spectrum is in the top trace. The real (middle 
trace) and the imaginary spectra (bottom trace) are obtained 
in a single shot AAQST experiment. The bar plots corre- 
spond to theoretically expected states (top row) and those 
obtained from AAQST experiments (bottom row). Fidelity 
of the AAQST state is 0.95. 



Ui n t(r) = exp(— iHt) is the unitary operator for evo- 
lution under the internal Hamiltonian H (see eqn. (1)) 
for a time r, and (^) are global x- rotations. The con- 
straint matrix M had 63 columns corresponding to the 
unknowns and 160 rows corresponding to the real and 
imaginary parts of 80 spectral lines. After optimizing 
the durations by minimizing the condition number us- 
ing a genetic algorithm, we obtained C(M) = 14.6 for 
n = 431.2 us and t 2 — 511.5 its. Again we study 
AAQST on two states: (i) Thermal equilibrium of the 
19 F spins: p\ = (a\ + a\ + uf)/2, and (ii) a ran- 
dom density matrix p 2 obtained by applying unitary 

Uq = (f ) x T o( 7r )^ T o (f )y % j witn T o = 2.5 ms, on ther- 
mal equilibrium state, i.e., p 2 = UopiU^. In both the 
cases, we initialize the ancilla i.e., the 1 H qubits on to a 
maximally mixed state by first applying a (ir/2) H pulse 
followed by a strong PFG in the z-direction. 

The real and imaginary parts of the single shot AAQST 
spectra, along with the reference spectra, are shown in 
Figs. 5 and 6 respectively. Again the line intensities 



?(i) 



{Rjjj v > i S)v\iv' } are obtained by curve-fitting, and 



all 



the 63 unknowns of the 3-qubit deviation density ma- 
trix are obtained by solving the matrix eqn. (5). The 



reconstructed density matrices along with the theoreti- 
cally expected states (pi and p?) are shown below the 
spectra in Figs. 5 and 6. The fidelities of experimen- 
tal states with the theoretically expected states (pi and 
P2) are respectively 0.98 and 0.95. The lower fidelity in 
the latter case is mainly due the imperfections in the 
preparation of the target state p2- The overall poorer 
performance in the liquid crystal system is due to the 
lower fidelities of the QST pulses, spatial and temporal 
variations of solute order-parameter, and stronger deco- 
herence rates compared to the isotropic case. In spite of 
these difficulties, the three-qubit density matrix with 63 
unknowns could be estimated quantitatively through a 
single NMR experiment. 



IV. CONCLUSIONS 

Quantum state tomography is an important part of 
experimental studies in quantum information process- 
ing. The standard method involves a large number of 
independent measurements to reconstruct a density ma- 
trix. The ancilla- assisted quantum state tomography in- 
troduced by Nieuwenhuizen and co-workers allows com- 
plete reconstruction of complex density matrix with fewer 
experiments by letting the unknown state of the input 
register to interact with an ancilla register initialized in 
a known state. Ancilla registers are essential in many 
of the quantum algorithms. Usually, at the end of the 



quantum algorithms, ancilla is brought to a state which 
is separable with the input register. The same ancilla 
register which is used for computation can be utilized for 
tomography after the computation. The ancilla register 
can be prepared into a maximally mixed state by dephas- 
ing all the coherences and equalizing the populations. 

We provided methods for explicit construction of to- 
mography matrices in large registers. We also discussed 
the optimization of tomography experiments based on 
minimization of the condition number of the constraint 
matrix. Finally, we demonstrated the experimental 
ancilla-assisted quantum state tomography in two sys- 
tems: (i) a system with two input qubits and one ancilla 
qubit in an isotropic medium and (ii) a system with three 
input qubits and two ancilla qubits in a partially oriented 
medium. In both the cases, we successfully reconstructed 
the target density matrices with a single quadrature de- 
tection of transverse magnetization. The methods intro- 
duced in this work should be useful for extending the 
range of quantum state tomography to larger registers. 
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